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SUMMARY 

The equilibrium equations and the compatibility conditions are fundamental 
to the analyses of structures. However, anyone who undertakes even a cursory 
generic study of the compatibility conditions can discover, with little effort, 
that historically this facet of structural mechanics had not been adequately 
researched by the profession. Now the compatibility conditions (CC's) have 
been researched and are understood to a great extent. For finite element dis- 
cretizations, the CC's are banded and can be divided into three distinct cate- 
gories: (1) the interface CC's, (2) the cluster or field CC's, and (3) the 

external CC's. The generation of CC's requires the separating of a local re- 
gion, then writing the deformation displacement relation (ddr) for the region, 
and finally, the eliminating of the displacements from the ddr. The procedure 
to generate all three types of CC's is presented and illustrated through exam- 
ples of finite element models. The uniqueness of the CC's thus generated is 
shown. 

The utilization of CC's has resulted in the novel integrated force method 
(IFM). The solution that is obtained by the IFM converges with a significantly 
fewer number of elements, compared to the stiffness and the hybrid methods. 

INTRODUCTION 

In the analyses of structures both the conditions of equilibrium and of 
compatibility come into play, except for the trivial statically determinate 
case. However the conditions of equilibrium are the most familiar to struc- 
tural analysts, perhaps because of a concern about the internal forces required 
for design and the wide acceptance of equilibrium as a universal and natural 
concept. In contrast, the concept of compatibility is much less familiar to 
structural analysts. The compatibility conditions (CCS's) were not known until 
mathematicians formulated them about a century ago (ref. 1), long after the 
equilibrium equations (EE's) had been derived. Even so, the early forms of the 
compatibility conditions were developed mainly for the manual analyses of sim- 
ple structures and were based on the concept of redundant structural elements. 
The notion of cutting the redundant members, which leads to the conditions of 
compatibility that have to be restored (henceforth referred to as classical 
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compatibility conditions (CCC's)), was formulated in the precomputer era. 
Formulating CCC's for large-scale computations proved inconvenient and ineffi- 
cient, so they have almost disappeared from current structural engineering 
practice. The general, and mathematically rigorous, compatibility conditions 
(analogous to St. Venant's strain formulation in elasticity (ref. 1)) have been 
formulated for finite element analysis. These finite element CC's are referred 
to as the global compatibility conditions (GCC's) or simply CC's. The GCC's 
are banded, and they are amenable to computer automation. (In stiffness analy- 
ses, linkage and continuity of displacements at element interface are popularly 
referred to as the compatibility conditions; in a strict sense, such con- 
straints are continuity conditions.) 

The utilization of compatibility conditions has resulted in the integrated 
force method (IFM) (refs. 2 to 16). In the IFM all the internal forces are 
treated as the primary unknown, and the system equilibrium equations are cou- 
pled to the global compatibility conditions in a fashion paralleling approaches 
in continuum mechanics, such as the Bel trami-Michell formulation in elasticity 
(ref. 17). 

The compatibility conditions, in basic form, have been introduced and com- 
pared with the CCC's in references 2 to 4. The purposes of this report are to 
(1) describe the physical aspects of the compatibility conditions, including 
the Interface CC's of finite element models; (2) demonstrate the generation of 
GCC's from their local counterparts; and (3) illustrate the benefits that 
accrue from the use of the global compatibility conditions In finite element 
analysis. 

The subject matter of this report is presented in the subsequent five sec- 
tions. In the second section the governing equations of the IFM are presented. 
In the third section we demonstrate the procedure to generate the compatibility 
conditions, and these concepts are illustrated in the fourth section. In the 
fifth section, comparison of results obtained by the IFM, the stiffness method, 
and the hybrid methods are presented, and the conclusions are given in the 
sixth section. 


EQUATIONS OF THE INTEGRATED FORCE METHOD 


In the Integrated force method, a discretized structure for the purpose of 
analysis Is designated as structure (n,m) where (1) structure denotes the types 
or structure (truss, Frame, plate, shell discretized by finite elements, their 
combinations, etc.) and (2) n and m are force and displacement degrees of 
freedom fof and dof, respectively. The structure (n,m) has m equilibrium 
equations and r = n - m compatibility conditions. The m EE's, C B ] { F } 

» {P}, and the r CC's, CC3CG3 {F} = {SR}, are coupled to obtain the governing 
equations of the IFM as follows: 


[B] I 
[C][G], 


{F} 



where 

[ B ] m x n equilibrium matrix 

[C] r x n compatibility matrix 


or CS]{F) = {P}* 


( 1 ) 
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[G] n x n concatenated flexibility matrix (containing material 
properties) and it links deformations {|3} to forces {F} as 
{13} = UGHF} 

{P> m-component load vector 

{6R} r-component effective Initial deformation vector such that 

{SR} = -[C]{|3q} where {|3q} Is the n-component Initial deformation 
vector 

[S] n x n IFM governing matrix 

The matrices [B], [C], [G], and C S ] are banded, and they have full row ranks of 
m, r, n, and n, respectively. 

The solution of equation (1) yields the n-forces {F}. The m-dl splacements 
{X} are obtained from the forces {F} by back substitution (ref. 9) as 

{X} = [J]{[G]{F} + {|3 0 }} (2) 

where [J] is the m x n deformation coefficient matrix defined as 
[J] = m rows of [[SH] T . 

Equations (1) and (2) represent the two key relations of the IFM for fi- 
nite element analysis that are needed to calculate forces and displacements, 
respectively. 


GENERATION OF THE COMPATIBILITY CONDITIONS 


The compatibility conditions and the associated coefficient matrix [C] are 
obtained from St. Venant's strain formulation In elasticity (ref. 1) as an ex- 
tension to discrete structural mechanics. The strain formulation Is illus- 
trated through the plane stress elasticity problem. The strain displacement 
relations (SDR) of the problem are 


e x “ 3x 

•y ' Ty > <3a) 

9u 3v 

Y xy * 3y + 8x^ 

Since in the SDR three strains (e x , ey, yxy> are expressed as functions of two 
displacements (u,v), the SDR contains one compatibility condition, which can be 
obtained by eliminating the two displacements from the three SDR as 


2 2 2 
a e a e a y xv 

^ + 8x 2 ~ ax 3y = 0 


(3b) 


The two steps that are necessary to generate CC's from St. Venant's strain 
formulation are (1) establish the strain displacement relations given by equa- 
tion (3a) and (2) eliminate displacements from the SDR to obtain the compati- 
bility condition given by equation (3b). 
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In the mechanics of discrete structures, the equivalent of SDR are the de- 
formation displacement relations (DDR). (Deformations of discrete analysis {(3} 
are analogous to strains {e} of elasticity.) The DDR can be assembled directly 
or obtained on an energy basis by utilizing the well known equality relation of 
internal strain energy and external work, which can be written in the case of a 
discrete structure (n,m) as 


\ {F} T {(3} = \ {X} T {P} (4a) 

Equation (4a) can be rewritten by eliminating loads {P} in favor of internal 
forces {F} by using the EE ([B]{F} = {P}) to obtain 

\ {X} T [B]{F} = \ (F } T {(3} or \ {F} T {[B] T {X} - {3}} = 0 (4b) 

Since the force vector {F} is not a null set, we finally obtain the following 
relation between member deformations and nodal displacements: 

{3} = CB] T {X} (5) 

The expression given by equation (5) represents the global deformation dis- 
placement relation applicable to any finite element model whose system equilib- 

rium equations can be symbolized as CB]{F} = {P}. In the DDR, n-deformations 
{3} are expressed in terms of m-dl splacements {X}; thus, there are r = n - m 
constraints on deformations, which represent the r compatibility conditions 
of the structure (n,m). The r CC's, in terms of total deformations, can be 
obtained by the elimination of m-di splacements from the n DDR, and in matrix 
notation the CC's can be written as 

CC]{3} = {0} (6a) 

The CC's, in terms of elastic deformations {3} that are given by equation 
(6b), are obtained from equation (6a) and fronrthe definition of the total de- 
formation {3}, which is composed of the initial deformations {3} 0 and the 
elastic deformations {3} 0 as {3} = {3} 0 + {3}g. Thus u 

[C] {3} g = {SR} (6b) 


where 

{SR} = -CC]{3 0 } 

[C] r x n global compatibility matrix 

The efficient generation of the CC matrix [C] is the subject matter of this 
report. The matrix [C] is rectangular and banded, and it has full row rank r. 
The CC's are kinematic relationships that are independent of sizing design 
parameters (such as area of bars, moments of inertia of beams, etc.), material 
properties, and external loads. The compatibility conditions depend on the 
initial deformation in the structure. For numerical efficiency, directly 
eliminating displacements from the DDR to obtain CC's is not recommended for 
large-scale computations. Instead, the global compatibility matrix [Cl is 
efficiently generated by utilizing such physical features of the compatibility 
conditions of finite element models as bandwldths, the determinacy of the grid 
points, and so on. 


4 


PROCEDURE TO GENERATE GLOBALCOMPATIBILITY CONDITIONS 

To generate the CC's, (1) separate a local region from the structural 
model on the basis of Interface, cluster, or external bandwidth considerations 
as explained later In this section; (2) establish the local deformation dis- 
placement relations (ddr) for the local region, and eliminate the displacements 
from the ddr to obtain the CC's for the region under consideration; and (3) re- 
peat steps (1) and (2) until all r CC's of the structure (n,m) are generated. 
These steps are elaborated in the section Illustrative Examples. The order of 
generation of the CC's is Immaterial; however, we recommend generating the in- 
terface CC's first, since these are most numerous, followed by the cluster 
CC's, and finally, the external CC's. 

Bandwldths of the Compatibility Conditions 

The CC's of discretized structures are banded. On the basis of bandwidth 
considerations, the CC's are divided into three distinct categories: (1) in- 

terface compatibility conditions, (2) cluster or field compatibility condi- 
tions, and (3) external compatibility conditions. By assuming the example of 
a finite element model as shown in figure 1, the three types of compatibility 
conditions can be illustrated as in figure 2. 

Interface Compatibility Conditions 

Numerous Interfaces internal to the structure are created in the discre- 
tization processes. The Interface is the common boundary shared by two or 
more elements. In the model shown In figure 1, the common boundary along 
nodes 1 and 7 is the interface between elements 1 and 2, the boundary connect- 
ing nodes 12 and 17 is the Interface between elements 13 and 14, and so on. 

The Interface between elements 1 and 2 is shown in figure 2(a). The deforma- 
tions of elements 1 and 2 must be compatible along the common boundary defined 
by nodes 1 and 7, which gives rise to Interface compatibility conditions. The 
number of CC's at the Interface depends on the element types (such as membrane, 
flexure or solid tetrahedron, etc.) and element numbers. The maximum bandwidth 
of the Interface compatibility condition can be calculated as 

JT 

H B «,CC-I>V <7a> 

J-l 


where 

MBWi cc maximum bandwidth of the Interface compatibility conditions 

JT total number of elements present at the Interface 

fof e j force degrees of freedom of the element j present at the Interface 

The bandwidth MBWj CC represents the maximum bandwidths of the Interface 
compatibility conditions written either in terms of forces { F} , as in 
[C][G]{F) = {0} (here we are referring to the bandwidth of the product matrix 
[C1CG], or in terms of deformations {|3}, as in [C]{|3} = {0} (here the 
bandwidth is that of the compatibility matrix [C]). The actual bandwidth of 
the compatibility matrix [C] is smaller than its maximum bandwidth. The 
Interface compatibility conditions of discrete analysis are analogous to the 
boundary compatibility conditions in elasticity (ref. 8). The interface CC's 
are the most numerous 
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in any finite element model. These are generated by writing the deformation 
displacement relation for the local region (such as shown In fig. 2(a) for the 
interface defined by nodes 1 and 7) and then eliminating the displacements from 
the local ddr as explained In the section Illustrative Examples. For the in- 
terface shown In figure 2(a), there are two elements (i.e., JT = 2). Let us 
assume that both are membrane elements; the fof of the triangular element 
foft Is 3 and that of the quadrilateral element, fofq Is 5. Then the MBW^ CC 
calculated from equation (7a) is 8. 

Cluster Compatibility Conditions 

Consider any element in the model shown in figure 1, such as 19. Ele- 
ment 19 along with its eight neighboring elements are shown in figure 2(b). 

The deformations of element 19 must be compatible with its neighboring elements 
(14 to 16 and 18 to 23). The CC's of the cluster of elements are referred to 
as the cluster compatibility conditions, which essentially represent the field 
CC's of St. Venant's strain formulation. The maximum bandwidth of the cluster 
CC's can be calculated as 


MBW 


ccc 


JTC 

E 

J-i 


(fof ej ) 


(7b) 


where 

MBW CC c maximum bandwidth of the cluster compatibility conditions 
JTC total number of elements present in the cluster 

Assuming that quadrilateral elements have fofq = 5 and triangular ele- 
ments have foft = 3, the MBW CC of the cluster shown in figure 2(b) can be 
calculated from equation (7b) as 41. 

In a finite element idealization, the number of cluster compatibility 
conditions are fewer than the number of Interface CC's. Generating cluster 
compatibility conditions requires the establishment of the ddr for the local 
cluster (such as shown In fig. 2(b) for element 19) and then elimination of 
the displacements from the local ddr as explained in the section Illustrative 
Examples. 


External Compati bi 1 i ty Condi tions 

Reactions are Induced at the nodes where displacements are restrained. If 
such restraints are sufficient only for the kinematic stability of the struc- 
ture, then the structure is externally determinate, and it has no external com- 
patibility conditions. If, however, the restraints on the boundary exceed the 
number of rigid body motions of the structure, then the structure is externally 
Indeterminate. The degree of external indeterminacy can be calculated as 
fol lows : 


R ext - N x - N f 


(7c) 


where 


6 



R ext number of external indeterminacy 

N x number of displacement components suppressed on the boundary 

Nf number of boundary conditions required only for the kinematic stability 

of the structure 

Let us assume that the finite element model shown in figure 1 represents a 
membrane structure. Then its external indeterminacy R ex t =7-3=4, since 
the number of actual boundary restraints N x is 7 and the kinematic stability 
requirement Nf is 3. 

To calculate the bandwidth of the external compatibility conditions, sepa- 
rate the local region connecting any two boundary nodes. Let the number of 
elements between the two nodes be represented by JTE, then the bandwidth of the 
external CC's is given by 


MBW 


ecc 


JTE 

£ 


(f0 V 


<7d) 


where 

MBW ecc maximum bandwidth of the external compatibility conditions 

Assuming as before that the quadrilateral elements have fof q = 5 and the 
triangular elements have foff = 3, the bandwidth of the external CC's shown in 
figure 2(c) can be calculated from equation <7d) as 8. 

If the boundary represents a determinate boundary, then no boundary CC's 
will be generated. The boundary CC's are obtained by eliminating the displace- 
ments from the deformation displacement relations written for the local bound- 
ary segment (e.g., for the model shown in fig. 1, the segment containing nodes 
1 and 3 and elements 2 and 3, also shown in fig. 2(c)) by following the proce- 
dure explained in Illustrative Examples. 

The Interface, cluster, and external CC's represent the local CC's. All 
three categories of local CC's are concatenated to form the system or the glo- 
bal CC's of the structure (n,m). The sum r = r^ cc + r ccc + r ecc of the num- 
ber of Interface CC's, cluster CC's, and external CC's is equal to r = n - m 
of the model. The values of n cc , r ccc , and r ecc can be calculated for 
discrete models; however it is not necessary to determine their values before 
generating the CC's. 


ILLUSTRATIVE EXAMPLES 

Examples of a few structures that are idealized by triangular membrane 
elements and bar elements are presented to illustrate the generation of global 
compatibility conditions from the local conditions such as interface CC's, 
cluster CC's, and external CC's. In the examples, triangular elements given 
by Przemieni ecki (ref. 18) and standard bar elements, which are adequate to 
illustrate the compatibility concepts yet simple enough for closed-form 
presentation, are chosen. The elements are shown in figure 3. The membrane 
element has three force unknowns, F] e , and ^3e! its s ^ x displacement 

degrees of freedom are X] e , X2 e . . . . , Xg e (fig. 3(a)). The bar element 
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has one force and four displacement unknowns (fig. 3(b)). For the membrane 
element, the 6 x 3 equilibrium matrix [B] e and its symmetrical flexibility 
matrix [G] e of dimension 3x3 are obtained in closed form (ref. 18) as 


-&12 

0.0 

*3l' 

~ m 12 

0.0 

m 31 

a i2 

-!l 23 

0.0 

m 12 

" m 23 

0.0 

0.0 

a 23 

- a 31 

0.0 

m 23 

~ m 31. 


(8a) 


where ijj and m^j denote the direction cosines of the direction vector de- 
fined by the edge ij; and 



sin 0 


3 


sin ©^ sin ©2 


cos 0 2 cot ©2 - v sin ©2 


cos © 2 cot ©2 - v sin ©2 


sin © 


sin © 2 s i n © 3 


cos ©^ot ©j - v sin © 1 


cos © 3 cot ©^ - v sin ©^ 


sin ©, 


cos © 2 cot ©2 - v sin ©2 cos © 3 cot © 3 - v sin © 3 


sin © 3 s 1 n © 1 


(8b) 


where 


E Young's modulus 

v Poisson's ratio 

t membrane thickness 


The angles ©i are defined in figure 3(a). 

Even though the 6x6 elemental stiffness matrix for the membrane can be 
generated in closed form, its explicit form is too complicated for presentation 
here. 


Example I - The Membrane 

Generation of the interface compatibility conditions Is illustrated 
through the example of a membrane shown in figure 4(a). The membrane is made 
of steel, with Young's modulus E = 30 ksi and Poisson's ratio v = 0.3; 
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dimensions a and b are 100 in. and thickness t = 1 in. The membrane is 
subjected to concentrated loads. The example Is also solved herein by the IFM, 
and the stiffness solution In symbolic form Is Included for comparison. The 
membrane is idealized by two triangular elements. The discretization has six 
force degrees of freedom, which are the concatenation of the two element force 
unknowns such that 


{F} T =■ <F 1 


lei 




where the subscript iej indicates the i th force of element j. 


(9a) 


A corresponding deformation component 3k. is associated with each force 
component F|<. The six-component deformation vector of the membrane is 


{3} = (3-| - 3 1 g -j i 1^2 = ^2el * • ■ ■ * 3g — 3^2^ (9b) 

The five system displacement degrees of freedom shown in figure 4(a) are rep- 
resented by {X} as 


{X} T = <X ] , X 2 X 5 > (10) 

The membrane is designated "membrane (6,5)" since it has six force and five 
displacement unknowns. Membrane (6,5) has five EE's and one CC ( 1 . e . , 
r = (n - m) . (6-5) = 1). 


Equilibrium equations of the membrane (6,5) . - The five system EE's of 
membrane (6,5), in terms of the six forces, are obtained from the elemental 
equilibrium matrices (refer to eq. (8a)) by following finite element assembly 
technique as 


'1 .0 

0.0 

0.7071 

0.7071 

0.0 

0.0 

-0.7071 

-0.7071 

0.0 

1 .0 

0.7071 

0.7071 

0.0 

0.0 

0.0 

0.0 

.0.0 

0.0 

0.0 

0.0 




r n 




F. 

1 

0.0 

0.0‘ 


1 


^50.00 

-1 .0 

0.0 


F 2 


100.00 

0.0 

0.0 


F 3 

= 

50.00 

• 1 .0 

0.0 


F 4 

r > 

50.00 

0.0 

1 .0. 


F 5 


100.00 





■s. J 


Since the system equilibrium matrix [B] has the dimension 5x6, the five EE's 
(eq. (ID) cannot be solved for the six forces; one CC is required to augment 
the five EE's to a solvable set of six equations in six unknowns. 


Compatibility condition for membrane (6,5) or the Interface compatibility 
condi tions . - Membrane (6,5) has one CC; therefore its local CC and global CC 
are represented by the same equation. The first step In the generation of the 
CC is to establish the deformation displacement relation for the membrane. The 
6x5 global DDR is obtained from equations (5) and (11) as 
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r > 

®i 


' 1.0 

0.0 

0.0 

0.0 

0.0' 


V 

*2 


0.0 

0.0 

1 .0 

0.0 

0.0 


X 2 



0.707 

-0.7071 

0.7071 

0.0 

0.0 



< 3 

> - 






< 

X 3 

*4 


0.707 

-0.7071 

0.7071 

0.0 

0.0 




0.0 

-1.0 

0.0 

1 .0 

0.0 


X 4 

w 


. 0.0 

0.0 

0.0 

0.0 

1.0. 


X 5 

< J 


Even though the single CC for this simple problem can be obtained by 
direct elimination of the five displacements from the six ddr given by equa- 
tion (12), this procedure Is not recommended because it can become numerically 
expensive for large structures. The concept of node determinacy, which greatly 
enhances computational efficiency in the generation of CC's, is presented next. 

The Node Determinacy Condition 

Node determinacy for general application Is presented first; then it Is 
specialized for membrane (6,5). Forces in determinate structures can be ob- 
tained from EE's alone; such determinate forces do not participate In the CC's. 
The concept of determinate structure Is extended to the nodes or grid points 
of a finite element jnodel , and to enhance computational efficiency, determinate 
nodes or grid points = are identified and eliminated at intermediate stages of 
the generation of CCs. 

Take any node of a finite element model, for example the node 1. Let 
K} represent the number of force components present in the EE's written for 

node 1. Let L} represent the displacement degrees of freedom of the node 1, 

which also Is the number of EE's that can be written at that node. The inde- 

terminacy of node i, designated (NR}), is defined as 

NR i = K 1 - L. (13) 

If NR} » 0, then node i is designated determinate. Forces present at a de- 
terminate node 1, referred to as determinate forces, do not participate In the 
compatibility conditions since such forces can be determined from the nodal 
equilibrium equations alone. Consequently, for determinate node 1, K} forces 
along with L} EE ' s T which correspond to L} = K} displacements, can be 
dropped simultaneously from the equilibrium matrix CB] to obtain the reduced 
equilibrium matrix [B]^ r ^ without affecting the CC's in any manner. Dropping 
of forces and displacements Is also equivalent to the elimination of appro- 
priate columns and rows of the deformation displacement relations. The reduced 
deformation displacement relation (designated DDR r }) that is obtained after Im- 
posing the node determinacy condition has the following form: 

DDR r] => {|3} (r1) = [B] (rl)T {X} <r1) (14) 

In equation (14), matrix [B] (r1)T has a dimension of {(n - K}) • (m - K})}. 

The deformation vector has the dimension n - K} , and displacements 
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{X}< r l> are of dimension m - . As expected, the number of CC's contained In 

the DDR r ] given by equation (14) Is r = n - m or r = {(n - K} ) - (m - Kj}, 
since no CC has yet been generated. The node determinacy condition has reduced 
the number of deformation displacement relations from n to n - K \ ; however, 
the number of compatibility conditions remains the same. The motivation behind 
dropping determinate variables at the Intermediate stage of the generation of 
the CC Is to enhance node determinacy at as many grid points as possible. 


For the example of membrane (6,5), observe that the last two of its EE's, 
given by equation (11), have two determinate forces, F 5 and Fg, and that 
these correspond to the two displacements X4 and X5 at node 4; therefore 
K]=2 and Li = 2. The node 1=4 Is determinate since NRi = Kj - Li = 0. 
The reduced DDR r i that Is obtained by taking into consideration the determinacy 
for node 4 has the following explicit form: 


r ^ 

*1 


' 1.0 

0.0 

0.0 ' 

3 2 

L _ 

0.0 

0.0 

1 .0 

3 3 

r 

0.707 

-0.707 

0.707 

P 4 


.0.707 

-0.707 

0.707. 

L J 






(15) 


The number of CC's given by the DDR r i (eq. (15)) still remains one (r = m - n 
= 1), since no CC has yet been generated. The local structure that Is obtained 
after the elimination of the determinate node 4 Is shown in figure 4(b). Since 
node 1 is fixed and node 4 has been dropped, the local ddr r ] given by equa- 
tion (15) corresponds to the DDR of nodes 2 and 3, which represents the Inter- 
face between elements 1 and 2. The Interface DDR has four deformations (3], 

32, 33 , and 34) expressed in terms of three displacements (X] , X2, and X3). 
The elimination of the three displacements from the four DDR given by equa- 
tion (15) yields the only CC associated with the Interface of the two membrane 
elements 1 and 2 along their common boundary defined by nodes 2 and 3). It has 
the following explicit form: 


(3 3 - iy = 0 (16) 

The CC given by equation (16) represents the deformation balance condition 
along the Interface of adjoining elements. Such CC's are referred to as the 
Interface CC's. For the membrane model, each Interface has one Interface CC. 
However the number of CC's at an Interface of any discrete model will depend 
on the type and number of elements that are connected to the interface. For 
example, the flexure problem given In the section Benefits Derived from the 
Compatibility Conditions has three CC's at each Interface. The Interface CC's 
represent the majority, though not all of the CC's of a finite element Ideal- 
ization. Generating Interface CC's requires that the ddr for the local inter- 
face be established and then that the displacements be eliminated from it. 

The CC, CC] { 3 } - {0}, for membrane (6,5), which Is given by equation (16), 
has to be expressed In forces, CC3CG] {F} = {0}, so that it can be coupled to 
the EE's (eq. (ID) to obtain the IFM governing equations, CS] {F} = {P}*. De- 
formations { 3 } are transformed Into forces {F} by the constitutive relation 
{3} = EG] {F} . Here, the 6x6 matrix EG] Is the block diagonal concatenation 
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of element matrices {G e } (see eq. ( 8 b)) for elements 1 and 2. The Interface CC 
(&3 - 34 = 0 ) In terms of forces has the following form: 


C 

F 1 


[0.5 0.5 1.0 


■ 1 .0 


-0.5 


-0.53/ 


> = (0} 


(17) 


For membrane (6,5) the maximum bandwidth of the CC (MBWj cc ) Is 6. The actual 
bandwidth of the compatibility matrix CC 3 (BW ac tual C* 2- The actual CC 
bandwidth of the composite matrix EC 3 [G 3 (BW ac tual CG^ 

The Integrated force method solution for membrane (6,5) . - For 
completeness, the solution of membrane (6,5) by the IFH Is presented. The IFM 
governing equation CS 3 {F} = {P} * Is obtained by coupling the EE's (eq. (ID) 
to the CC (eq. (17)) as 
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The solution to equation (18) yields the six forces. The five displacements 
are obtained from the forces by back substitution in equation ( 2 ) as 
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Analysis of the membrane by the stiffness method . - For comparison, the 
membrane Is also analyzed by the stiffness method. The stiffness equations are 
well known but complicated; therefore, the analysis is carried out in symbolic 
form. To establish the paralleism between the IFM and the stiffness method, a 
slightly different procedure than normal is followed; the purpose will become 
evident in the process of the solution. For the membrane, a displacement vec- 
tor {X c } of dimension 12 that represents the concatenation of the 2 elemental 
displacement degrees of freedom (see fig. 4 (c)) is defined as: 

{ V = < X cl = X lel’ X c 2 = X 2 e 1 ’ ' ' • ’ X c 1 2 = X 6 e 2 ) (20) 

where the subscript iej represents the 1 th displacement component for the 
j th element. 

Notice the similarities between the displacement vector {X c } of equa- 
tion ( 20 ) and the force vector { F } given by equation ( 9 a). These vectors rep- 
resent the concatenation of the elemental displacements and the elemental force 
degrees of freedom, respectively. The equilibrium matrix in the stiffness 
method, In terms of concatenated displacement vector {X c }, is obtained by fol- 
lowing assembly techniques as 

[CK 1 ] : CK 2 ]]{X c } = {P} (21) 

The stiffness matrix CK-j ] has the dimension 3x6. Its three rows represent 
contributions to the EE's at node 2 along X], and at node 3 along X 2 and X 3 
(see fig. 4(a)). Likewise, the matrix [K 2 I has the dimension 5x6, which con- 
tributes to the equilibrium at node 2 along X], node 3 along X 2 and X 3 , and 
node 4 along X 4 and X 5 . Equation (21), which represents 5 EE's in terms of 
12 unknown displacements cannot be solved for the 12 variables. Seven dis- 
placement continuity conditions are required to augment the EE's given by 
equation (21) to arrive at a solvable 12 x 12 system. The seven displacement 
continuity conditions of the membrane are as follows: 


X 1 e 1 

= 0 

X 2 el 

= 0 

( 3el 

= X 3e2 

: 4el 

X 4e2 

: 5el 

X le 2 

X 6 e 1 

= 0 

X 2 e 2 

= 0 i 

J 


( 22 a) 


The seven displacement continuity conditions can be represented by a single 
matrix equation 
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[CTY]}X} = {0} (22b) 

where the matrix CCTY] Is the displacement continuity matrix of the dimension 
7 x 12. 

The displacement continuity conditions that are given by equation (22) and 
the EE's in terms of displacements that are given by equation (21) are coupled 
to obtain the 12 x 12 equation system (eq. (23)) from which the 12 concatenated 
displacements X c can be obtained: 


■[K ] :K 2 I 

«C> - < 

r *\ 

P 



. CCTY] . 


0.0 

k J 


From the displacements, the internal forces can be calculated by back 
substitutions . 

From the structures of IFM equations (eqs. (1) and (18)) and from the 
equation of the stiffness method (eq. (23)), we observe the following: ' (1) In 
the IFM, the EE's, written in terms of forces, are augmented by the CC's, also 
in terms of forces; (2) in the stiffness method the EE's are expressed in terms 
of displacements and these EE's are augmented by the displacement continuity 
conditions; and (3) the IFM equations (eq. (18)) are fewer in number, and also 
sparser than the stiffness equations (eq. (23)). Details about equation spar- 
sity and the computations required to generate the solution by the IFM and the 
stiffness method are given in reference 14. The IFM satisfies both the EE's 
and the CC's simultaneously, whereas the stiffness method is based on the EE's 
and displacement continuity conditions. 

Example II - The Two-Bay Membrane 

The generation of both the interface and the cluster CC's are illustrated 
by using the example of a two-bay membrane shown in figure 5. The membrane is 
discretized by eight triangular elements, and it has eight nodes. The model, 
designated "membrane (24,15)," has n = 24 force unknowns and m = 15 dis- 
placement unknowns. Membrane (24,15) has m = 15 EE's and r = n - m = 9 
CC's. Because of the increase in complexity, the algebra for the example is 
presented in symbolic form. The system equilibrium matrix [B m ] of dimension 
15 x 24 for membrane (24,15) is assembled by following standard techniques. 

The global deformation displacement relations (GDDR) that correspond to the 
equilibrium matrix [B m ] can be symbolized as 

GDDR => {B} - CB m ] T {X} (24) 

The GDDR (given by eq. (24)) contains 9 CC's, since its 24 deformations are 
expressed in terms of 15 displacements. The model shown in figure 5 has eight 
interelement boundaries, defined by nodes 1 and 5, and nodes 5 and 8, and so 
on. From example I, we know that each interface has one CC; therefore, the 
eight Interfaces yield eight CC's, which can be generated by following both 
parts of step I, explained in the following paragraphs. Steps 1(a), 1(b), II, 
and III, are parts of the general procedure to generate the CC's. For clar- 
ity, the steps are explained by using the example of membrane (24,15) as an 
1 1 lustration. 
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Step 1(a) - Local structure and Interface compatibility conditions . - Con- 
sider any Interface - for example, the Interface along nodes 1 and 5 (between 
elements 1 and 2) for membrane (24,15). Separate the Interface and the ele- 
ments as shown In figure 6 (a). The local structure shown therein is statically 
unstable; therefore Impose adequate numbers of restraints to make the local 
structure kinematically stable. The restraints do not influence the CC's, and 
they can be Imposed at any of the nodes of the local structure. The structure 
requires 2 restraints, which are imposed at node 5. The stable local struc- 
ture (6,5) has one interface CC, designated as CC^. The deformation displace- 
ment relation for the local structure, designated ddr |_^ , is extracted from the 
global GDDR. The local ddrjj consists of six deformations (3], 32. • • • . 3fi) 
expressed in terms of five displacements (X] , X 4 , X 5 , Xg, and X 7 ). The single 
compatibility condition in ddru , is generated by following the procedure given 
for example I, and it turns out to be 33 - 34 = 0. 

Step 1(b) - Update the global deformation displacement relations . - The 
number of GDDR of the structure (n,m) is reduced after the generation of each 
CC. The reduced GDDR, which is designated GDDR r ] , has mi = m rows and 
n l = n - n C ] columns, where n C ] represents the number of CC's generated in 
step 1(a). For the example in step 1(a), n c i is 1 since only one CC was 
generated In Step 1(a). The row and column dimensions of GDDR r i are 
ni = (24 - 1) = 23 rows and mi = 15 columns. The GDDR r i contains eight CC's 
since only one out of r = 9 CC's of structure (24,15) has been generated. 

The GDDR r i Is obtained by dropping one deformation displacement relation from 
the GDDR. Any deformation that has participated in the CC generated in step 
1(a) can be dropped. For the example, deformation 33 or 34 can be dropped. 

Step II - Local structure and its interface compatibility conditions . - 
The local structure consisting of elements 2 and 3 is separated next, and its 
Interface CC's along nodes 2 and 6 are generated by following steps 1(a) and 
1(b). Steps I and II are repeated until all the Interface CC's have been gen- 
erated. For this problem there are eight Interface CC's at the end of whose 
generation the GDDR r s will have 16 rows and 15 columns containing 1 CC. 

Steps I and II are sufficient to generate all the Interface CC's, which 
are the most numerous CC's in all finite element models. Since the interface 
CC depends on the few elements that are common to the Interface, the computa- 
tion time required for the generation of such CC's after the equilibrium matrix 
has been established is insignificant compared to the total solution time. 

Step III - Cluster compatibility conditions . - In a finite element model, 
a cluster Is defined as a series of adjoining elements. The cluster compati- 
bility conditions represent constraints on the deformations of the elements 
that belong to the cluster. A cluster can be generated for any element. For 
membrane (24,15), take element (4). The cluster for this element, shown in 
figure 6 (b), consists of four elements, 2 to 4 and 8 , and six nodes, 1 to 3, 

5, 6 , and 9. Let us designate its deformation displacement relation as ddrQj . 
The cluster is stable, so there is no need to impose any of the restraints in- 
dicated In step 1(a). If the cluster was unstable, then it would be necessary 
to Impose the additional constraints indicated in step 1(a). The ddrcu of the 
cluster contains nine deformations, 34 , 36. to 3 i 2 . ancl 322 to 324 
(note that 35 , 37 , and 3 s have been eliminated during the generation of the 
Interface CC), which are expressed in terms of nine displacements, X] to X 3 , 
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Xg to Xg, X^ 4 , and X 15 . The cluster, which has nine deformations expressed 
In nine displacements Is determinate and contains no additional CC's. 

The second cluster that Is defined for element 2 Is shown In figure 6 (c). 
The ddr^ of the second cluster contains 7 elements and 14 deformations (since 
7 deformations have been dropped during the generations of the 7 Interface CC's 
belonging to the cluster), expressed In terms of 13 displacements. The cluster 
contains one CC, which Is obtained by following step 1(a): 

{0.447(|3 6 - 0 2 + |3g - 3 20 > - 0 ]6 - 0 ]2 + 0.894(0, - P 5 + 0 ]? - 0^)} = 0 


(25) 


Step 1(b) should be exercised to reduce the GDDRd 2 - The GDDRqi 2 has 13 de- 
formations expressed In terms of 13 displacements, and It contains no CC's, 
thereby Indicating that all 9 CC's (8 Interface CC's and 1 cluster CC) have 
been generated for membrane (24,15). 

Example 3 - The Stiffened Membrane 

The generation of external CC's and Interface CC's when different element 
types are present In the discretization Is Illustrated by taking the example of 
the stiffened membrane shown in figure 7(a). This membrane Is discretized by 
8 membrane elements and 16 bar elements. Nodes 1 and 3 of the membrane are 
fully restrained. The fof of the membrane (n = 40) consists of the 24 mem- 
brane forces and 16 bar forces. The dof of the membrane Is m = 14. The 
membrane Is designated "membrane (40,14);" It has 14 EE's and r = (40 - 14) 

= 36 CC's. The global deformation displacement relation for the structure, 
designated as GDDR sm and consisting of 40 deformations expressed In terms of 
14 displacements, Is assembled by following standard techniques. 


Interface compatibility conditions of membrane (40,14) . - The first Inter- 
face, which Is defined by nodes 1 and 5, and associated membrane elements 1 
and 2, and bars 9 to 12 and 16, is considered. One boundary constraint, the 
displacement at node 2 , Is Imposed for its overall stability as shown In fig- 
ure 7(b). This local structure has njj = 11 fof consisting of six membrane 
forces and five bar forces, and mu =5 dof. The local structure, shown in 
figure 7(b), Is designated as SMLjj (11,5); Its six (rjj = 6 ) Interface compat- 
ibility conditions, which can be generated by following step I, are: 


(0 3 
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P 2 - 
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Pl 6 = ° 


* 10 -°} 


&12 - 0 
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(26) 


( 27 ) 
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The CC given by equation (26) represents the membrane-to-membrane Inter- 
face compatibility condition, and those given by equation (27) represent the 
five membrane-to-bar Interface CC's. 

Steps I and II are repeated for Interfaces 4 and 5, 4 and 8 , 5 and 8 , 

5 and 9, 5 and 6 , 2 and 6 , and 2 and 5 to generate, respectively, 3, 3, 3, 3, 

2 , 3, and 1 additional local Interface CC's. After the generation of the 
24 Interface compatibility conditions, the reduced deformation displacement 
relation has 2 CC's and consists of 16 deformations that are expressed In 
terms of 14 displacements. 

Cluster and external compatibility conditions of membrane (40,14) . - Step 
III Is evoked and yields two compatibility conditions: one Is a cluster CC 

Identical to the CCC of membrane (24,15) as given by equation (25), and the 
other 34 + 3 ]g = 0 , which represents a constraint on the deformations between 
boundary nodes 1 to 3 , Is the external compatibility condition of the membrane. 
Generation of the external CC is further explained with the example of a bridge 
truss. 


External compatibility conditions of a bridge truss . - A bridge truss sup- 
ported at two nodes that are far apart (see fig. 8 (a)) illustrates the genera- 
tion of external CC's. The bridge truss being analyzed can be designated as 
truss (26,20); It has 20 EE's and 6 CC's, and its global GDOR has 26 
deformations expressed In terms of 20 displacements. Skeletal structures such 
as trusses and frames do not have any Interface CC's; their CC's can be either 
cluster or external ones. 


Cluster compatibility conditions of the truss (26,20 ). - The cluster for 
element 1, consisting of six bars, is shown In figure 8 (b). The unstable 
cluster Is made stable by Imposing a constraint at node 3. The cluster, desig- 
nated "bay (6,5)," has one CC, which is obtained by following step I. At this 
time any one bar, for example bar 1 that corresponds to deformation 3]. is 
dropped; the resulting bay (5,5) is shown In figure 8 (c). Steps I and II are 
repeated until all five cluster CC's are generated and the structure shown in 
figure 8 (d) is obtained. 


External compatibility conditions . - The reduced structure shown in fig- 
ure 8(d) has 1 CC since Its ddr consists of 21 bar deformations expressed In 
terms of 20 displacements. The single CC is obtained by first imposing anode 
determinacy condition that reduces the ddr to six deformations expressed in 
five displacements and then by eliminating the displacements. The external CC 
(32 + Bio + Bi 5 + 320 + 325 = represents a homogeneous constraint on all 
the deformations between boundary nodes 1 and 12, as shown In figure 8(e). The 
bandwldths of the external CC's for restrained nodes that are far apart, typi- 
cally encountered in long span bridges, can be large from physical considera- 
tions, since the deformations between fully constrained node of a truss that 


are far apart have to be zero 




The larger bandwldths of few 


external CC's do not Impose any major problem because In the IFM the solution 
process Is carried out by using sparse matrix techniques. Quite often such ex- 
ternal CC's can be trivially generated by mere Inspection. 


17 


Uniqueness of Compatibility Matrix 

The compatibility conditions are homogeneous equations; therefore the 
CC's can be multiplied by any nonsingular matrix, for example CR U ] , to obtain 
the feasible compatibility matrix CCu ] given by equation (28). The feasible 
matrix CC U 3 is a linear combination of the rows of matrix CC]: 

CR U ][C]{3} * { 0 } or [C u ]{ 3 > = { 0 } (28) 

The procedure presented in this paper generates the matrix CC] and not CC U ], 
This can be proved by observation. Take the example of the two membrane inter- 
face CC's (CC] and CC 2 ) defined as: 


CC 1 

* P kl 

- P k 2 " 0 

(29a) 

CC 2 

* p k3 

- 3 k4 = 0 

(29b) 


A linear combination of CC] and CC 2 yields a feasible compatibility condition 
CC 3 as 


CC 3 * B kl - B k 2 * S k3 - »k4 - 0 (29c) 

Notice that in equation (29c), deformations of CC] such as 3^ and 3^2 
and deformations of CC 2 such as 3|<3 and 3k4 are present. However, after 
the generation of CC] , one of the two deformations 3^1 or 3^2 must be 
dropped (step 1(b)); therefore their combination cannot occur in subsequent 
CC's. In other words, after CC] has been generated, the feasible compatibility 
condition CC 3 , which includes CC] cannot be obtained. Dropping a deformation 
that has participated in the CC, immediately after its generation, avoids the 
possibility of its further participation in any other CC. The process gener- 
ates the matrix CC ] and not its combinations such as the usable CC matrix C C u ] . 
Therefore the unique C C ] matrix is generated. 

Product of Compatibility and Equilibrium Matrices - a Null Matrix 

The product of the CC matrix CC] and [EE] matrix [B] Isa null matrix 
(CC3CB3* = [0]). This can be verified by direct substitution of the DDR in the 
CC's as 


CC3 {3} = {0} (30a) 

and 

{ 3 } = CB3 T {X} (30b) 

Next, eliminate deformations in favor of displacements between equations (30a) 
and (30b) to obtain 


cb] t cc] = cc: t cb] = [0] ooc) 

The null product property of the two fundamental EE and CC operators implies 
that errors in equilibrium equations can propagate to the compatibility con- 
ditions and vice versa. 
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BENEFITS DERIVED FROM COMPATIBILITY CONDITIONS 


The accuracy of solutions obtained by the Integrated force method, the 
stiffness method, or the hybrid method Is of paramount Importance, since all 
are approximate formulations. All methods attempt to satisfy the equilibrium 
equations written In terms of forces or displacements. However, explicit com- 
patibility conditions, In a strict sense, are Imposed only In the Integrated 
force method. Any Improved solution accuracy, of the Integrated force method 
over the other formulations should be a consequence of the explicit presence 
of the global compatibility conditions in the IFM. Based on the theory of the 
integrated force method, a Generalized Integrated Force Technique (GIFT) com- 
puter code has been developed. To illustrate the solution accuracy in finite 
element calculations, a plate bending problem Is solved by the GIFT code and 
other standard analyzers such as MSC/NASTRAN (ref. 19), ASKA stiffness codes 
(ref. 20), a mixed formulation MHOST (ref. 21), and Chang's hybrid method 
(ref. 22). 

The plate parameters considered (see fig. 9) are the following: size of 

the plate; a = b = 40 In. (101.6 cm); aspect ratio a/b, varied between 1 and 2; 
thickness of plate, h = 0.2 in. (5.08 mm); Young's modulus, E = 30 000 ksi 
(21 091.81 kg/mm^); Poisson's ratio, v = 0.3; and magnitude of concentrated 
load at the center, P = 500 lb (226.795 kg). Both simply supported and clamped 
boundary conditions are considered. 

To compare solution accuracy, the problem Is solved by using two types of 
elements: a four-node rectangular element and a three-node triangular element. 

The IFM elements assume three forces (such as a shear force and two bending 
moments) and three displacements (a transverse displacement and two slopes) per 
node, as depicted in figure 9. A cubic polynomial with 12 constants is used to 
approximate the transverse displacement In the element field. Normal moments 
M x and My are assumed to have linear distributions, and the twisting moment 
M X y is constant In the element domain (ref. 14). 

The elements of the general purpose programs NASTRAN, ASKA, and MHOST are 
specialized to generate only the flexure solution. The elements used are the 
following: (1) QUAD4 (both ASKA and MSC/NASTRAN have QUAD4 elements), 

(2) TRIB3 (triangular element of ASKA program), (3) TRIA3 (triangular element 
of MSC/NASTRAN), and (4) TUBA3 of ASKA code (a higher order triangular element 
with six dof per grid point). For the first three elements, which are well 
known in the literature and popular In practice, the bending response of the 
elements represents three dof per node. 

The hybrid elements have more unknowns, for example, for flexural re- 
sponse. Chang's program has the equivalent of seven unknowns at the nodes, 
whereas the mixed formulation MHOST has more unknowns per grid point. The IFM 
elements and the stiffness elements (such as QUAD4 of MSC/NASTRAN, QUAD4 of 
ASKA, TRIB3 of ASKA, and TRIA3 of MSC/NASTRAN) are equivalent with respect to 
their nodal degrees of freedom. The hybrid and TUBA3 elements are higher order 
elements than those of the IFM. 

In the stiffness method, nodal stress parameters that are calculated from 
displacements by back substitution are discontinuous and ambiguous at grid 
points (ref. 23); therefore calculation of forces at the nodes are routinely 
avoided In the stiffness method. In this situation the noncontroversial nodal 
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displacement is used in the comparison. Remember however that in the IFM, 
forces are the primary variables from which the secondary displacement un- 
knowns are obtained by back substitution, The central deflection of the plate 
w c , given by Timoshenko (ref. 24), is 0.2036 in. (5.715 mm). 

MacNeal (ref. 25) Introduced a grading scheme for the evaluation of finite 
elements as fol lows : 

A <2-percent error 

B 2- to 10-percent error 

C 10- to 20-percent error 

D 20- to 50-percent error 

F >50-percent error 

The results obtained by all four formulations were graded by using MacNeal 's 

scheme and are presented in tables I to III. Results obtained by the IFM and 
the stiffness formulations are also presented graphically In figures 10 and 11. 
The IFM results for simply supported and clamped boundary conditions for dif- 
ferent aspect ratios are presented in table IV. From the numerical results of 
the plate flexure problem presented in tables I to IV and figures 10 and 1 1 , we 
observe the following: 

(1) For the IFM rectangular element, convergence occurs for the first 
model, which consists of four elements. If symmetry Is taken into considera- 
tion, then convergence occurs for a single-element idealization. Both stiff- 
ness (MSC/NASTRAN and ASKA) and hybrid methods (MH0ST and Chang's) converge 
slowly. To achieve an A grade, MSC/NASTRAN QUAD4 element idealization requires 
36 elements, whereas ASKA QUAD4 secures only a B grade, even for 100 elements. 
The hybrid method of Chang secures an A grade for 64 elements, whereas MHOST 
secures a B grade for the same level of discretization. 

(2) For the IFM triangular elements (see fig. 11), the result Is dis- 
cernible from an analytical solution for the first model, which has four ele- 
ments, but even so, the result displays engigeerlng accuracy. The next model, 
with eight elements, converges to the analytical solution and also achieves an 
A grade. None of the stiffness elements, such as TRIA3 of MSC/NASTRAN and 
TRIB3 and TUBA3 of ASKA, could secure a grade of A, even for models with fine 
discretization. 

(3) The IFM result for a simply supported boundary follows the pattern of 

a clamped boundary; namely. It secures a grade of A for the first model, which 

has four elements. The IFM element retains an A grade for aspect ratios up to 

1.6, but for the ratio 2.0, a total of eight elements, which corresponds to a 

2x4 mesh, Is required to secure an A grade. Other examples more or less fol- 
low the pattern depicted in tables I to IV. Overall the IFM convergence rate 
Is very fast whereas both the stiffness and hybrid methods converge slowly or 
struggle to do so. 


CONCLUSIONS 

1. The structural mechanics profession recognizes that both equilibrium 
equations (EE's) and compatibility conditions (CC's) are essential for stress 
analysis. However, the compatibility conditions In typical finite element cal- 
culations were promoted via such concepts as cutting and closing the gaps, or 
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displacement matching at nodes (deflection or slope should have unique values) 
and so on. Although such concepts are somewhat related to the CC's of finite 
element models, they do not represent the true CC's that are analogous to the 
strain formulation of St. Venant. The true compatibility conditions of finite 
element analysis have been understood to a great extent, though we do not claim 
that the understanding Is total. Attempts should be made by the profession to 
understand the CC's in totality, rather than to avoid them because they are 
mathematically formidable or analytically more difficult than the familiar 
equl 1 i bri urn equations . 

2. In finite element analysis, the system equilibrium equations In terms 
of forces or displacements can be assembled from element matrices. The ques- 
tion Is, can such an assembly technique be developed for the compatibility con- 
ditions also? The generation of compatibility conditions is not equivalent to 
the direct-assembly technique of the finite element analysis, even though there 
is a close resemblance in that the global compatibility conditions are assem- 
bled from their local counterparts such as Interface CC's, cluster CC's, and 
external CC's. We do not yet know a direct assembly scheme, but such a possi- 
bility has not altogether been ruled out. We do, however, believe that the 
compatibility generation scheme given in this paper is rather elegant, since 
element characteristics, connectivities, and such, which are already contained 
in the equilibrium equations (and consequently in the deformation displacement 
relation because it is the raw Ingredient of CC's) are referred to only once. 

3. The compatibility conditions are unique. The computation time required 
to generate the CC's is a small fraction of the total solution time. 

4. The quality of the solution in approximate methods is dependent on the 
extent to which equilibrium equations and compatibility conditions are satis- 
fied. Since the integrated force method (IFM) satisfies both the EE's and the 
CC's simultaneously, the solution via the IFM is accurate, as expected. The 
stress parameters obtained by the stiffness and hybrid methods do not satisfy 
the EE's even at the cardinal grid points; therefore solution quality by such 
methods Is prone to be poor in comparison to IFM results. 

5. Since all the finite element formulations are approximate in nature, 
we recommend generating solutions via the integrated force method, and by the 
stiffness method and then comparing them, rather than qualifing the results 
obtained by any one formulation by successive mesh refinements. 
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APPENDIX - SYMBOLS 
dimensions of rectangular element 
equilibrium matrix of dimension m x n 
element of equilibrium matrix 
system equilibrium matrix 
compatibility matrix of dimension r x n 
compatibility condition 
classical compati bi 1 i ty condi tion 
displacement continuity matrix 
local deformation displacement relation 
displacement degrees of freedom 
Young's modulus 
equi 1 i bri urn equations 
force vector of dimension n 
force degrees of freedom 
flexibility matrix of dimension n x n 
element flexibility matrix 
global compatibility conditions 
global deformation displacement relations 
Integrated force method 

deformation coefficient matrix of dimension m x n 

total number of elements at interface 

total number of elements in cluster 

total number of elements between two boundary nodes 

stiffness matrix of dimension m x m 

direction cosines along edge ij 

maximum bandwidth 

plate bending moments 
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m 

NRj 

n 

{P} 

r 

[S] 

SDR 

t 


u,v 


w c 

{X} 

{X c } 

{3} 

3k 

Yxy» E x> c y 


number of displacement variables 
Indeterminancy of node 1 
number of force variables 
load vector of dimension m 
number of compatibility condition; (n - m) 
I FM governing matrix of dimension n x n 
strain displacement relations 
membrane thickness 
membrane displacement components 
central deflection of plate 
displacement vector of dimension m 
concatenated displacement vector 
deformation vector of dimension n 
k^h deformation component 
strain components 
Poisson's ratio 


Subscripts : 

ccc cluster compatibility conditions 

ej element j 

ecc.ext external compati bi 1 i ty condl tlons 

1 cc interface compatibility conditions 

q quadrilateral 

t triangular 

0 initial 


Superscripts : 

r number of compatibility conditions 

rl reduced 

T transpose 
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TABLE I. - REPORT CARD FOR RECTANGULAR ELEMENT 
WITH CLAMPED BOUNDARY CONDITION 


Number of elements 
for full plate (mesh) 

IFM 

MSC/NASTRAN 

QUAD4 

ASKA 

QUAD 4 

4 (2x2) 

A 

F 

f 

: 

16 (4x4) 

A 

B 

B 

36 (6x6) 

A 

A 



64 (8x8) 

— 

A 



100 (10x10) 

- 

- 

> 
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TABLE II. - REPORT CARD FOR INTEGRATED FORCE 
METHOD WITH RECTANGULAR ELEMENT 


Number of elements 
for ful 1 pi ate (mesh) 

Aspect 

ratio 

Clamped 

boundary 

Simply 

supported 

boundary 

4 (2x2) 

1 .00 

A 

A 

16 (4x4) 

1 .00 

A 

A 

4 (2x2) 

1 .20 

A 

- 

4 (2x2) 

1 .40 

A 

- 

4 (2x2) 

1.60 

A 

- 

4 (2x2) 

1 .80 

B 

- 

4 (2x2) 

2.00 

B 

- 

8 (2x4) 

2.00 

A 

- 


TABLE III. - REPORT CARD FOR RECTANGULAR ELEMENT 


Number of elements 
for full pi act (mesh) 

IFM 

Mi xed 

method 

MHOST 

Hybrid 
method 
Chang 
(ref. 22) 

4 (2x2) 

A 

F 

F 

16 (4x4) 

A 

C 

C 

64 (8x8) 

A 

B 

A 


TABLE IV. - REPORT CARD FOR TRIANGULAR ELEMENT 


Number of elements 
for ful 1 plate 

IFM 

MSC/NASTRAN 

TRIA3 

ASKA 

TRIB3 

ASKA 

TUBA3 

4 

B 

F 

_ 

F 

8 

A 

D 

F 

F 

16 

A 

C 

- 

- 

32 

- 

B 

C 

D 

128 

- 

- 

B 

B 
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(a) Interface CC. 
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Figure 2. - Bandwidth of compatibility conditions (CC’s). 
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(a) Membrane element; F^, F 2> and F3 are force degrees of freedom {fof) and X^ t X 2 , . . Xg 
are displacement degrees of freedom (dof). 


Y 


X 4 



Fl 


0 

(b) Bar element; F 1 is fof and X-j, X 2 , X^, and X4 are dof. 
Figure 3 . - Membrane and bar elements. 
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(a) Membrane. 


(b) Interface nodes 2 and 3. 
Figure 4. - Analysis of membrane. 






(c) Conatenated model. 


Figure 5. - Two -bay membrane. 



(a) Interface CC. (b) Determinate duster. 

Figure 6. - Compatibility conditions for two-bay membrane. 



(c) Indeterminate cluster. 



(b) Interface CC. 


Figure 7. - Composite (membrane and bar) structure. 
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